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13  ABSTRACT 


The  problem  of  determining  the  optimal  inspection  epoch  is  studied  for  relia¬ 
bility  systems  in  whic.«  N  components  operate  in  parallel.  Life  time  dis¬ 
tribution  is  arbitrary  but  known.  The  optimization  is  carried  with  respect 
to  two  cost  factors:  the  cost  of  inspecting  a  component  and  the  cost  of  failure. 
The  inspection  epochs  are  determined  so  that  the  expected  cost  of  the  whole 
system  per  Lime  unit  per  cycle  will  be  minimized.  The  optimization  process 
depends  in  the  general  case  on  the  whole  failure  history  of  the  system.  This 
dependence  is  characterized.  The  cases  of  Weibull  life  time  distributions 
are  elaborated  and  illustrated  numerically.  The  characteristics  of  the  optimal 
inspection  intervals  are  studied  theoretically. 
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1. _ Introduction. 
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In  the  present  study  we  investigate  the  problem  of  determining 
the  optimal  inspection  epochs  of  a  reliability  system  which  is  comprised 
of  N  components,  operating  independently  (in  parallel)  and  having  the 
same  life  time  dis  buttons.  The  life  time  distribution  is  known. 

An  inspector  visi*'t  the  system  at  a  predetermined  inspection  epoch  and 
finds  a  certain  n*  fiber  of  components  which  have  failed.  The  exact  times 
of  failure  are  "  ■-aown.  All  the  components  v;hich  have  failed  during 
the  interval  be  jp.en  inspections  are  replaced  by  new  components.  Com¬ 
ponents  which  ,j  a  not  failed  are  left  in  the  system.  We  consider  two 
types  of  cost  'actors:  (i)  The  cost  of  inspection,  which  depends  on 
the  number  ol  components  in  the  system;  and  (ii)  The  cost  of  failure 
per  unit  tiro.1.  This  cost  component  measures  the  loss  due  to  a  failure 
of  components.  The  objective  is  to  determine  an  inspection  policy  that 
would  be  optimal  with  respect  to  the  criterion  of  minimising  the  total 
expected  (discounted)  cost  for  the  entire  future.  However,  since  we 
are  <.ieal:.;;g  with  cases  of  general  life  time  distributions  'not  neces¬ 
sarily  exponential)  the  dynamic  programming  solution  is  excessively 
complicated,  even  in  the  truncated  case  (when  the  number  of  inspections 
should  not  exceed  a  prescribed  bound).  Therefore,  we  are  considering 
in  the  present  paper  a  sequential  myopic  procedure.  Accordingly >  after 

each  inspection  the  epoch  of  the  next  inspection  is  determined,  as  a 
function  of  the  whole  past  failure  history  of  the  system.  The  aim  is 
to  minimize  the  conditional  expected  cost  per  time  unit  from  the  pre¬ 
sent  time  until  the  next  inspection  epoch.  In  the  case  of  exponen¬ 
tial  life  time  distributions  (constant  failure  rates)  the  optimal 
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inspection  interval  -(time  interval  between  inspections)  does  not  depend 
on  the  past  history  of  the  system.  As  shown  in  the  present  study  if 
the  life  time  distribution  is  not  exponential  this  dependence  might  be 
very  strong,  especially  if  N  is  not  large  and  the  life  time  distribution 
is  of  a  decreasing  failure  rate  (DFR).  The  dependence  of  the  optimal 
inspection  intervals  on  the  observed  number  of  failures,  and  on  the  num¬ 
ber  of  components  that  were  replaced  at  previous  inspections  and  are  still 
operating,  will  be  explicitly  characterized.  We  start  in  Section  2  by 
formulating  the  model  and  the  associated  distributions.  In  Section  3 
we  develop  a  general  formula  for  the  sequential  determination  of  the 
optimal  length  of  the  inspection  intervals.  In  Section  4  we  derive  the 
corresponding  formulas  for  life  time  distributions  of  the  Weibull  family; 
and  illustrate  the  process  with  a  numerical  example.  In  Section  5  we 
try  to  explain  the  complex  process  illustrated  in  the  example  of  Section 
4  by  further  theoretical  development. 

There  are  numerous  papers  in  the  reliability  literature  on  inspec¬ 
tions  epochs  and  optimal  maintenance.  For  the  general  theory  see  Chapter 
4  of  Barlow  and  Proschan  [1],  Articles  which  are  close  to  the  present 
study  are  those  of  Kamins  [4],  Kander  [J>],  Kander  and  Naor  [6]  and  Kander 
and  Rabinovitch  [7].  The  present  study  provides  further  elaboration  of 
a  chapter  in  the  thesis  of  Fenske  [3].  The  main  difference  between  the 
present  study  and  the  acticles  mentioned  above  is  in  the  basic  model. 

The  present  study  is  concerned  "•>  uh  mu I “icomponent  systems  while  che 
other  studies  treat  the  whole  -vs*- «.  as  one  component.  The  study  of 
Ehrenfeld  [2]  was  based  on  a  modal  similar  to  ours,  but  Ehrenfsld  considered 
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the  problem  of  determining  the  inspection  interval  for  the  estimation 
the  mean  time  between  failures  in  the  exponential  case. 


2.  The  model  and  associated  distributions. 

Consider  a  reliability  system  which  consists  of  N,  N  >  1,  com¬ 
ponents.  These  components  operate  independently  (in  parallel).  Let  T 
designate  the  life  time  of  a  component.  This  is  a  random  variable  having 
a  known  distribution  function  (c.d.fl  F  (t).  We  assume  that  F(t)  is 
absolutely  continuous,  with  a  positive  density  function  f(t),  0  <  f(t)  <  0, 
and  F(0)  =  0.  We  further  assume  that  the  expected  value  of  T,  according 
to  F(t)  is  finite.  Let  SQ  s  0  and  let  S0  <  Si  <  Sg  <  • • •  <  ^  <  • • • 
designate  a  sequence  of  inspection  epochs.  Let  Jm  (m  =  1,2,...)  designate 
the  number  of  components  that  failed  during  the  time  interval  (S  , ,  S  ). 
All  the  Jm  components  are  replaced  at  the  inspection  epoch  S^.  The 

N  -  J  components  which  have  not  failed  during  (s  ,  >  S  )  are  classi¬ 
cs  m-1  m 


(m)  .(m) 


(m) 


fied  into  m  disjoint  subsets  A  Q,  A  ...,  A^_^.  The  subset  A^ 

(j  =  0,  m-1)  contains  all  the  components  that  were  replaced  at  epoch 

S.  and  did  not  fail  throughout  the  time  interval  [S.,  S  ).  Let  n^ 
j  J '  m  J 

designate  the  number  of  elements  of  A^.  Obviously,  A^m|^  c  A^ 

and  <  n^  for  each  j  =0,  1,  ...  and  m  =  j,  j  +  1,  ...  Let 

(m)  r  j  (m)  /  (m)  (m)  (m).  , 

n  =  J  ,  and  nN  '  -  (n  n .  n  ,  ,  ...,  n'  )  for  eacn  m  =  0,  1,  ...: 

mm  ~  01m 

»(0>E  »• 

If  a  component  belongs  to  the  subset  A^  then  its  conditional 
life  time  distribution  at  time  t  is: 

,(m) 


(2.1) 


F'7  (t)  =  P(T  s;  t-s ,|t  ?  S  -S.l 
J  1  J  m  y 


,  if  t  s  s 


m 


F(t-s,)-F(S,-s|)|  lf  t  >  s 


m 
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In  particular,  F^m\t)  =  F(t-S  ),  if  t  >  S  and  zero  otherwise, 
m  m  ra 

The  conditional  densities  of  U  =  T  -  (S  -S.),  corresponding  to  the  life 

®  J 

time  T  of  a  component  which  belongs  to  AV  play  an  important  role 
in  our  procedure.  We  can  call  U  the  remaining  life  time.  If  a  com¬ 


ponent  is  chosen  at  random  ft  time  t  =  S  +0  its  remaining  life  time 

m 

U  has  a  conditional  density  function 

(2  2^  h  fuls^  _  I  T'  n<m)  f(U+Sm  V  „  >  0 

(2.2)  hm(ul£  )  -  N  n  j  l-F(Sm-Sj)  ' 


where  =  (S,,  ...,  S  ). 

~  i  m 

We  notice  that  if  T  has  «'  negative  exponential  distribution  i.e., 
f(t)  =  \e  **t,  t  a  0,  for  any  0  <  \  <  co,  then  h  (u|s^m\  n^)  =  f(u) 
for  all  m  =  1,  2,  ...  and  all  (S^m\  n^).  This  is  a  well  known 
property  of  the  negative  exponential  distributions.  Let  H  (ul  S^.n^) 
designate  the  c.d.f,  corresponding  to  (2.2). 


.2: _ Sequential  determination  of  inspection  epochs. 

We  will  consider  in  the  present  section  the  problem  of  deriving  an 
inspection  policy  which  attains  a  certain  economic  objective.  We  assume 
therefore  that  the  cost  of  inspecting  the  system  is  $  CQ  per  inspection 
and  on  the  other  hand,  if  an  element  fails,  then  the  cost  associated  with; 
its  failure  is  $  Cf  per  time  unit.  The  inspection  policy  adopted  here 
is  the  following.  Given  the  history  connected  with  the  past  in  inspection 
intervals,  i.e.,  (S^  \  n^)  determine  the  (m+1)  st  inspection  epoch 

so  that  the  average  expected  cost  per  time  unit  of  inspection  and  of  fail¬ 
ure,  over  the  (m+1)  st  inspection  interval  will  be  minimized.  We  re¬ 
mark  in  this  connection  that  this  policy  is  in  essence  a  myopic  policy, 


which  minimizes  the  expected  time  average  costs  for  each  inspection 
interval  individually.  A  Dynamic  Programming  determination  of  the  in¬ 
spection  epochs  could  attain  a  more  global  optimization.  However,  attempts 
at  Dynamic  Programming  solutions  lead  to  complicated  sets  of  recursive 
functional  equations.  The  solution  of  these  equations  is  generally  very 
tedious. 


Let  A  designate  the  length  of  the  (m+1)  st  inspection  interval. 
That  is,  A  =  S  , ,  ~  S  .  Given  (S^m\  n^)  the  conditional  expected 

HIV!  IT 

average  cost  per  time  unit,  under  A>  is 


.v  p  # a  Q(m)  n(mK  _J3  C*  (m)  .  v  f^ufSm  Sj^ 

(3.i)  Rn  S  ,  n  >  =  a"  T  ^  nj  (A"U)  l-F(Sm“S^) 


Or  in  terms  of  the  conditional  distribution  of  the  remaining  life  time 
U  we  can  express  (3.1)  in  the  form 


(3.2) 


R  (AJs(m), 

rn 


n(m))  = 


r  +  NCfHm(A|S(m),  n(m)) 


~  N  f  uh  (u I s* 
A  •>  m' 


The  optimal  (m+1)  st  inspection  epoch  is  defined  as  S  .  =  S  +  A°, 

m+i  m 

where  A°  is  a  positive  real  value,  A,  for  which  the  infimum  of  (3.2) 
is  attained. 

Let 

(3.3)  p  “  J  uh  (u|S^m\  n^)  du, 

0 

be  the  expected  remaining  life,  given  (S^m\  n^).  According  to  the 

rsj 

assumption  of  the  vious  section,  <  «>.  Differentiating  R  (AjS^m\n^m^) 

with  respect  to  A  obtain  that  if  u  £  C  /NC,.  then  A°  ?  », 

mol 
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This  is  a  case  in  which  no  more  inspections  are  warranted.  On  the  other 
hand,  if  >  C  /NC^,  t*iere  exists  a  unique  solution,  A0,  to  the 
equation: 

(3.4)  f  uh  (uls^,  n^)  du  =  C  /NC,. 

M  in  1  rv*  0  1 

0 

We  realize  from  (3.4)  that  S  is  a  function  of  the  statistic  (S^m\n^m^) 
of  the  system. 

As  we  have  already  mentioned  in  cases  of  exponential  life  time 
distributions  the  optimal  length  of  the  inspection  intervals  is  the  same 
for  all  m  =  1,  2,  ...  .  If  0  =  X  *  is  the  mean  time  between  failures 
(MTBF )  in  the  exponential  case  then  =  9  for  all  m,  and  the  condi¬ 
tion  for  a  finite  A°  is  thatCQ  <  NCf9;  i.e.,  the  cost  of  inspecting 
an  element  is  smaller  than  the  expected  cost  of  failure  of  an  element. 

If  this  condition  is  satisfied  then,  letting  y  =  Cq/NC^9,  it  is  easy 
to  show  that 

(3.5)  A°  «  |  *  [4J  , 

V 

2 

where  [4]  designate  the  y-fractile  of  a  chi-square  distribution 

with  4  degrees  of  freedom. 


4,  Optimal  inspection  epochs  for  Weibull  distributions. 

Suppose  that  the  life  time  of  an  element,  T,  follows  a  Weibull 
distribution,  with  a  density  function 


(4.1) 


f (t;9,a) 


0  ,  if  t  s  0 

a  ta  1exp{-t0!/ 9}  ,  if  t  >  0; 
9 


where  a  and  0  are  positive  real  parameters.  We  notice  that  if  0  <  a  <  1 


then  the  distribution  has  a  decreasing  failure  rate  (DFR),  and  if  1  <  a  <  oo 
its  failure  rate  is  increasing  (IFR).  When  a  =  1  tl-e  distribution,  is 
exponential.  Given  (S^m\  n^)  the  density  function  of  the  remaining 
life  U  assumes  the  special  form 


(4.2) 


exp\(S  -S.)  /9>. 
L  m  3  l 


for  0  i  u  s  oo,  When  m  =  0  (4.2)  reduces  to  (4.1).  Following  the 
procedure  given  in  the  previous  section  we  realize  that  <  °°  if  and 
only  if, 

(4.3)  ro/N  <  cf  e1/a  r<£  +  D 

^1/a  +  the  expected  life  time.  If  (4.3)  is  satisfied  then 

the  optimal  value  of  is 

(4.4)  si  =  [|G_1(‘vi|,  D/1/0t  , 

where  G  \vjp,v),  is  the  y-fractile  of  the  Gamma  distribution  G(p,v), 
with  scale  parameter  p,  and  where  y  =  CQ/ (NC^B^^rX"  +  1)).  We  notice 
that  if  2/a  is  a  positive  integer  j  then 


We  determine  now  a  general  expression  for  the  left  hand  side  of 


(3.4).  According  to  (4.2), 


<M>  |  u.h<9’“>  (u|s<“>,  du  = 


i  ni  .  »  ^  . 

h  njm  exP^sm"sj)a/^-  e  I  u(u+Vsj)a-lexp{-(u+Vs.)^e3 


By  a  proper  change  of  variable  we  obtain 

«,7)  f  JU(U  +  S^S^axpHu^-S  )“/«}  du  . 
0 

<VV«a/e 

J  je1^  w1'  ^  -  (.>m-Sj)"j  exp{-w)  dw  = 


(S  -S.)a'9 
m  ] 


e1/aT(£  +  1}  •' 


as  -s.+a) 

m  i 


a 


L.  9 
I -  (S  o  )& 

(Sm-S.)  ;exp{ - ^  ■  '  }  -  exp{ 

L 


1  \ 
h  5  +  l) 


\  f<S„-S  >°  !  \1 

r1- ;  l*  s  -  ljj 


v.S  _S  +A) 

— 5B ] _ I 

0  j 


Substituting  (4.7)  into  (4.6)  we  cbtain  that  S  .  =  S  +  A.  where  A 

m+1  ra 

is  the  root  of  the  equation: 
ra 


sC* 


<M>  J  S  n(m)  oxPf(s  -s  )“/e)  ,d.  , 

j=o  j  31  9  L> 


!  (s  -S.+A)  .  \ 

' - ;  iA+i) 

“  J 


y  0  ,  '  a  y 

_ I _ I  r  n<m)  ^  o  x 

K1/ar(i)  Nj=0  J  “  J 


=  Y  + 


1  -  exp{- J  [(S  -  ?  '  '"a 


|  0  Lni  Sj+A) 


(S  -S.) 
m  j 


a 
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V  is  as  before. 

We  notice  that  for  m  =  0  the  solution  of  (4.8)  is  reduced  to 

the  one  given  by  (4.4).  In  Figure  1  we  illustrate  the  solution  of  (4.8) 

for  three  Weibull  distributions,  where  the  sequences  were  generated 

by  Monte  Carlo  simulation.  The  cases  under  consideration  have  the  following 

parameters:  Cf  =  $10.  ,  Cq  =  $200‘N  ,  0  =  [hr]  100  and  ft  =  3/4,  1 

and  5/4.  The  case  of  ft  --  1  corresponds  to  the  exponential  distribution 

with  mean  9  =  100,  According  to  (3*5)  the  optimal  inspection  interval 

2 

for  a  =  1  is  of  length  [hr.]  50  Xy  [4]  where  y  =  Co/N*Cf9  =  0.2  .  One 

, .  ,  o 

can  find  in  any  statistical  tables  that  [4]  =  I.65.  Hence,  the 

• 

optimal  interval  between  inspections  is  in  the  exponential  case  is  length 
82.5  hours.  The  case  of  ft  =  5/4  represents  an  IFR  distribution.  We 
see  in  Figure  1  that  the  optimal  inspection  intervals  are  of  length  which 
vary  very  little  around  59  hours.  It  is  interesting  to  notice  that  in 
the  present  case  of  an  IFR  distribution  the  optimal  inspection  intervals 
do  not  depend  strongly  on  the  number  of  components,  N,  in  the  system. 

This  is  not  the  case  when  the  Weibull  disti ibution  is  a  DFR  (ft=3/4). 

As  illustrated  in  Figure  1  the  optimal  intervals  for  DFR  distributions, 
as  obtained  from  (4.8),  are  sensitive  to  N.  When  N  =  10  there  are 
considerable  fluctuations  of  the  solution  of  (4.8).  When  N  =  100  these 
fluctuations  diminish.  The  general  trend  of  growth  in  the  length  of  the 
inspection  intervals  is  however  the  same.  An  explanation  of  thia  phenom¬ 
enon  will  be  provided  in  the  next  section.  Finally  we  remark  that  the 
numerical  solution  of  equation  (4.8)  in  the  case  discussed  here  has  been 
attained  followin'  the  Newton-Raphson  iterative  corrections  to  an  initial 
solution.  For  fu:.'.  •>  details  see  Fenske  [3]. 


Inspection  Interval  (hours) 
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!'’i  ure  1.  Option:”  Infection  Int- rvals  for  ..'uiiull 
Distribution  ./ ith  C,  «  <10,  CQ  =--  -20GK,  and  0  =  IGO(hrs) 
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The  variance  o£  W  is 


m 


f 


m 


\tn)(D(m))2 

(5.7)  Var[Wm3  =  ~  <  2 


\j=o  U-Hsm-sf¥ 


e(m)D(m) 

J  J 

1-F(S  -S .)/ 
m  jV 


2  \ 
J 


We  have  shown  that  foi  any  sequence  of  inspection  <  ochs  Var  {W^}  =  0(N*"^) 
as  N  -*  oo.  This  explains  why  the  fluctuations  oj.  the  roots  of  (4.8)  are 
relatively  larger  when  N=10  and  small  when  N=100.  We  consider  now  a 
particular  sequence  of  inspection  epochs  which  consists  of  values  of 
obtained  by  the  repeated  r  lution  (for  each  m)  of  the  equation  tnm  =  Y>  i.e., 


(5.8) 


u  ‘  I 

P  uf(u)du  +  Tj  1  - 

0  >l  L 


J2 


9<j) 


A 

j‘  uf  (u+S^-S  .)du 
0 


Y 


A 

is  the  root  A  of  ^  uf(u)du  =  y,  ■  and  for  each  m  =  1,  2,  . the 

(m+l)st  inspection  epoch  is  given  by  ^  +  A.  The  sequence  of 

fixed  inspection  epochs  determined  by  this  procedure  corresponds  to  the 

expected  values  of  n^  and  we  therefore  label  this  procedure  as  the 

Procedure  Of  Averages.  In  Table  1  we  provide  the  inspection  intervals 

determined  by  the  Procedure  of  Averages,  and  the  corresponding  multinomial 
/  \ 

probabilities  9^  '  (j=,  . ..,  m, ),  for  the  two  cases  represented  in 

Figure  1,  The  graph  of  the  corresponding  inspection  intervals  for  the 

case  of  a  =  3/4  (DFR)  is  also  plotted  in  Figure  1.  As  is  demonstrated 

in  Table  1,  in  the  IFR  case  (a=5/4)  the  significant  contribution  to  the 

solution  is  expected  to  be  that  of  n'  '  ana  n^  ,,  or  of  their  corre- 

m  m-1 

sponding  expected  values.  Furthermore,  the  optimal  length  of  the  inspec¬ 
tion  intervals  varies  very  little  with  the  number  of  inspections,  m,  and 
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1  -  F(Sm) 

j-1  <”'> 

(1  -  ZL  e,H;  )  [1  -  F(S  -s  )],  j  =  1,  ...,  m. 
i=u  i  m  J 

It  follows  that  for  any  fixed  sequence  of  inspection  epochs  and  for  each 
m  =  0,  1;  •  1 1 

£{n<m)}  =  N  9<m) 

L  J  J 

(5.3)  j  «  o,  1,  ...,  m 
Vai[a<B)}  =  N  9<B)<l-e<"°) 

and 

(5.4)  cov(n^m\  n^m^)  =  -  N0^  ofm\  all  0  <c  j  <  k  <.  m. 

J  K  J  K 


,(m) 


(5.2) 


and 


9i 


(m) 


From  (5.2)  and  (p.j)  we  conclude  that  if  the  length  of  each  inspection 

interval  is  not  smaller  than  A  then  for  any  distribution  F>  lim  n.  =  0 

o  J  m-fco  j 

for  each  j. 

The  variable  W  is  a  linear  combination  of  multinomial  random 
m 

variables.  Its  expectation  is 


(5.5) 


m 


to 


m 


E{w  }  =  D(m)  +  2 
'  m  o  ,  , 
J=l 


ur 


(m) 

D.  > 
J 


where 

(m)  _  n 

D .  u  u  f (u+S  -S  . )  du. 

J  0  m  j7 


(5.6) 
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its  expectation  reaches  in  the  present  example  a  stable  situation  after 
two  inspections.  This  is  not  the  case,  however,  in  the  DFR  distribution 
(a=3/4).  The  probabilities  9^  approach  zero,  as  m  grows,  '»ery  slowly. 
This  is  reflected  in  a  steady  increase  in  the  length  of  the  inspection 
intervals  as  m  grows,  and  a  stable  situation  is  reached  in  the  present 
example  only  after  10  inspections. 

To  insure  chat  the  inspection  intervals  discussed  in  Sections  3 
and  4  will  have  similar  properties  to  those  determined  by  procedures  of 
fixed  inspection  epochs  we  could  consider  the  following  adjustment. 

First,  determine  for  each  m  =  1,  2,  ...  two  fixed  sequences  of  inspec¬ 
tion  epochs  which  will  constitute  upper  and  lower  (confidence)  limits 
for  the  solution  of  (3.4)  (or  (4.8)).  This  can  be  done  by  utilizing 
formulae  (5-5)  and  (5«7).  The  lower  confidence  limits  could  be  obtained 
by  repeated  solution  (for  the  root  A)  of  the  equation 

(5.9)  U)m  +  3.[Var{Wra}]1/2  =  y  ;  m  =  1,  2,  ... 

The  upper  limit  can  be  obtained  by  solving  the  equation 

i  /p 

(5.10)  -  3.[Var{Wmj]  =y,  m  =  1,  2,  ... 

In  the  second  phase  of  computation  solve  equation  (3.4).  If  the 
solution  lies  between  the  roots  of  (5.9)  and  (5.10)  proceed;  otherwise 
truncate  the  solution  to  eithe  *  the  lower  limit  or  to  the  upper  limit, 
whichever  is  closer  to  the  actual  solution.  Such  an  adjustment  will 
guarantee  that  every  inspection  interval  will  be  bounded  by  lower  and 
upper  values  which  are  determined  by  fixed  sequences  of  inspection  epochs, 
and  will  therefore  have  general  characteristics  as  established  here. 


Table  1:  Values  of  optimal  inspection  intervals  A  [hrs]  and  multimonial 
probabilities  under  the  Procedure  Of  Averages  for  Weibull  distributions 
with  0  =  100  [hrs]  and  cost  :omponents  Cq  =  $200N,  =  $10. 


Case  IS  a  =  5/4  (IFR) 


2 

j 

=3 

j=4 

j=5 

j 

=6 

j=7 

j=3 

j 

=9  j=10 


0.2017  0.7983 


59.3  0.0211  0.1543  0.8246 


59.1  0.0016  0.0161  0.1600  0.8223 


0.0010  0.0013  O.OI67  0.1595  0.8225 


0  0.0001  0.0013  0.0166  0.1595  0.8225 


0  0.0001  0.0013  0.0166  0.1595  0.8225 


Case  ii:  a  =  3 A  (DFR) 


147.6  0.6548  0.3452 


156.8  0.4825  0.2216  0.2959 


164.8  0.2839  0.1231 


0.0311 


0.1880 

0.2830 

0.2771 

0.1039 

0.1302 

0.1742 

0.2735 

0.08  03 

0.0984 

0. 1267 

0. 1716 

0.0530 

0. 0761 

0. 0957 

0.1246 

0.0500 

0.0600 

0. 0739 

0.0941 

0.0401 

0. 0473 

0.0580 

0.0726 

fflSI 

0.0379 

o.ouc 

0.03  59 
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